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SECTION  I 


INTRODUCTION 


This  document  presents  a detailed  description  of  the  work 
performed  in  the  conduct  of  the  "Compressor  Cascades  with 
Splitter  Vanes  Analysis  Program"  under  USAF  Aerospace  Research 
Laboraties  (ARL)  Contract  No.  F33615-75-C-1004 . The  objectives 
of  the  program  were: 

o To  determine  the  most  desirable  configuration  of  a 
transonic  splitter  vane  in  a given  cascade 

o To  present  the  results  in  a manner  which  would  per- 
mit their  use  for  proper  splitter  vane  selection 
in  other  cascade  configurations 

These  objectives  were  pursued  by  the  inviscid  analysis  of 
transonic  flow  through  cascades.  Calculations  were  performed 
with  use  of  the  methods  defined  in  References  1 auJ  2.* 


The  Program  was  divided  into  the  following  four  phases: 

o Phase  I - Analysis  of  Cascades  Without  Splitter  Vanes 

o Phase  II  - Analysis  of  Cascades  With  Splitter  Vanes 

o Phase  III  - Cascade/Splitter  Vane  Design  Study 

o Phase  IV  - Generalization  and  Development  of  Design 

Rules 


The  addition  of  splitter  vanes  within  a cascade  is  used,  in 
general,  to  improve  the  turning  of  a blade  row.  However,  due 
to  the  added  surface  area  of  the  passage,  the  cascade  with 
splitter  vanes  may  be  expected  to  suffer  additional  viscosity- 
induced  losses.  At  present,  there  are  no  analytical  means  by 
which  the  designer  may  compare  the  potential  gain  in  turning 
with  the  higher  loss  of  the  cascade  so  that  the  best  splitter 
geometry  may  be  chosen  to  accomplish  a particular  performance 
goal.  When  the  complication  of  transonic  operation  is  added, 
the  designer  has  little  experimental  data  for  guidance. 

Significant  analytical  methods  used  in  each  phase  of  the 
study  are  discussed  in  Section  II,  while  the  actual  results  for 
each  phase  are  given  in  Section  III. 


‘References  follow  Section  IV. 
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SECTION  II 

TECHNICAL  DISCUSSION 


1.  PHASE  I - ANALYSIS  OF  CASCADES  V7ITHOUT  SPLITTER  VANES 

a.  Objectives 

The  objectives  for  Phase  I were  to  compare  analytically 
determined  blade  loadings  and  turning  angles  with  experimental 
data  for  the  cascade  tests  described  in  Reference  3.  The  basic 
analytical  tool  used  in  this  phase  was  the  inviscid  transonic 
flow  program  developed  by  AiResearch,  which  is  described  in 
Reference  1.  The  comparison  was  made  first  on  the  basis  of  the 
predictive  ability  of  the  numerical  method,  and  second,  on  the 
design  usefulness  of  the  results  of  such  a calculation. 

b.  Inviscid  Transonic  Flow  Calculation 

The  inviscid  transonic  flow  calculation  method  used  in  this 
investigation  required  that  transformation  of  the  potential 
equation  be  made  into  a variable,  non-orthogonal  grid  system 
within  the  cascade.  The  transformation  resulted  in  an  equation 
that  involved  both  elliptical  and  hyperbolic  operators.  The 
coefficients  of  each  of  the  operators  depend  upon  the  local  Mach 
number  and  the  local  grid  angle.  For  example,  when  the  flow  is 
supersonic  and  the  grid  angle  is  in  the  characteristic  direc- 
tion, the  coefficient  of  the  elliptic  operator  is  zero. 

The  extent  of  the  analysis  of  Phase  I w;js  limited  to  the 
following  cascade  flow  conditions: 

o Cascade  design  inlet  flow  conditions  (M.  = 1.46, 

= 66.85°  as  determined  by  experiment]. 

o Three  values  of  the  cascade  outlet  Mach  number, 

which  implied  three  different  values  of  the  overall 
static-pressure  rise. 


Specification  of  the  inlet  Mach  number,  inlet  flow  angle, 
and  outlet  Mach  number  for  an  inviscid  flow  implies  a relation 
between  the  stream  tube  contraction  across  the  cascade  and 
the  outlet  flow  angle.  This  relation,  as  shown  in  Equation  (1)  , 
is 


i 


1 

1 

i 
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where  the  b-width  is  defined  as  the  distance  normal  to  the  plane 
of  potential  flow  calculation  between  the  effective  flow  bounda- 
ries at  the  cascade  end  walls.  Equation  (1)  is  plotted  for  the 
specified  inlet  flow  conditions  in  Figure  1 for  a range  of  outlet 
Mach  number.  This  figure  shows  that  if  it  were  possible  to 
throttle  a cascade  through  the  subsonic  discharge  flow  range, 
and  maintain  a constant  outlet  flow  angle 

increases  as  M-  decreases.  (Note,  decreasing  M2  is  equivalent 
to  increasing  P2/Pj^.) 

The  b-width  ratio  can  be  related  by  the  continuity  equation 
directly  to  the  axial  velocity-density  ratio  (AVR)  across  the 
cascade . 


1 

AVR 


where 


AVR  = 


^1^X1 


The  experimental  observation  that  turning  increases  with  increas- 
ing AVR  (or,  alternately  expressed,  turning  increases  with 
decreasing  b2/bj^)  is  illustrated  in  Figure  2 (taken  from  Refer- 
ence 4)  for  a typical  low-speed  compressor  blade  element.  The 
same  trend  can  be  calculated  with  the  use  of  an  inviscid  flow 
analysis  coupled  with  an  assumption  concerning  the  Kutta  condi- 
tion. If  the  Kutta  condition  is  used  to  fix  the  circulation  of 
the  cascade,  reasonable  agreement  is  obtained  between  the  calcu- 
lations and  experiments  (Reference  4) . Two  different  inviscid 
calculations  are  illustrated  in  Figure  2 for  a cascade  having 
relatively  small  viscous  effects.  Experience  has  shown  that  when 
viscous  effects  become  more  important,  agreement  often  breaks 
down . 


At  least  two  ways  in  which  the  inviscid  transonic  analysis 
method  may  be  set  up  for  comparison  with  experimental  measure- 
ments may  be  used. 

(1)  The  uniform  flow  angle  ^he  far  downstream 

station  may  be  fixed  to  agree  with  the  experimental  value.  When 
doing  this,  the  normal  stagnation  point  would  change  from  that 
expected  from  the  Kutta  condition  at  the  trailing  edge.  To 
obtain  periodicity,  a stagnation  point  and,  therefore,  the  con- 
nection of  the  dividing  streamline  would  have  to  be  moved  to  the 
suction  surface.  This  method  was  eliminated  because  of  numer- 
ical difficulties  imposed  on  the  grid  system. 
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Figure  1.  Variation  of  the  Product  of  the  b-Width 
Ratio  and  the  Cosine  of  the  Outlet  Flow 
Angle  with  Outlet  Mach  Number. 
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1 Figure  2.  Turning  Angle  Versus  Meridional  Velocity 
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(2)  The  Kutta  condition  can  be  satisfied  at  the 
trailing  edge  by  suitable  choice  of  the  outlet  flow  angle.  Thus, 
an  inviscid  estimate  of  turning  is  obtained. 

The  second  method  is  more  suited  to  the  program  goals,  but 
the  first  method  often  leads  to  better  prediction  of  the  blade 
surface  pressure  distribution. 

c.  Cascade  Geometry  and  Test  Data 

A summary  of  pertinent  geometrical  details  of  the  cascade 
tested  in  Reference  3 is  presented  in  Table  I.  The  blade  ele- 
ment is  illustrated  in  Figure  3,  which  indicates  the  blade  sur- 
face coordinate  positions  chosen  to  numerically  represent  the 
blade  element  in  the  present  analysis. 

TABLE  I.  CASCADE  PHYSICAL  CHARACTERISTICS. 


Chord 

3.004  in. 

Axial  chord 

1.8397  in. 

Blade  spacing 

1.5810  in. 

Blade  span  - inlet  plane 

2.861  in. 

Blade  span  - exit  plane 

1.4"7  in. 

Maximum  thickness/chord  ratio 

0.03614 

Metal  angle  - leading-edge  pressure  surface 

61.417® 

Metal  angle  - leading-edge  suction  surface 

65.479® 

Mean  camber  angle  - leading  edge 

63.448" 

Mean  camber  angle  - trailing  edge 

22.534® 

Stagger  angle 

52.316® 

Camber  angle 

40.913® 

Solidity 

1.9 

The  end  walls  of  the  cascade  investigated  experimentally  in 
Reference  3 were  convergent,  so  that  the  overall  geometrical 
b-width  ratio  (bj/b.)  was  0.4871.  The  variation  of  the  geomet- 
rical b-width  with  axial  position  through  the  cascade  is 
illustrated  in  Figure  4. 

Typical  experimental  data  (from  Reference  3)  for  the  ARL 
supersonic  compressor  cascade  is  shown  in  Figures  5 and  6 in 
terms  of  wai’  ‘ — Ic-pressure  ratio  versus  axial  distance. 

Passage  No.  4 of  a seven-bladed  cascade  was  instrumented  with  10 
static  pressure  taps  on  each  sidewall  to  obtain  this  data. 
Pressure  ratio  is  expressed  as  wall  static  pressure  over  the 
upstream  total  pressure.  This  pressure  ratio  and  the  correspond- 
ing isentropic  flow  relation  were  employed  to  obtain  the  experi- 
mental wall  Mach  numbers  presented  in  Section  III. 
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AXIAL  LENGTH/AXIAL  CHORD,  X/C^ 


Figure  4.  Variation  of  the  Geometrical  b-Width  as 
a Function  of  Axial  Distance  Through  the 
Cascade. 


Figure  5.  Experimental  Data  for  Supersonic 
Compressor  Cascade  (Reference  3;. 
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d.  Design  Method 

An  inviscid  calculation,  no  matter  how  accurately  it  pre- 
dicts surface  pressure  distribution,  does  not  by  itself  directly 
provide  the  results  needed  for  design  decisions.  For  the  case  in 
point,  two  principal  aerodynamic  criteria  are  used  in  design. 

(1)  A measure  of  the  loss  (or  relative  position  on  a 
loss  scale)  of  the  section  relative  to  other  possible  sections. 

(2)  A measure  of  the  increase  in  turning  obtainable  by 
the  section. 

A useful  estimate  of  the  relative  performance  of  one  configura- 
tion versus  others  has  been  obtained  from  inviscid  calculations. 
Total  losses  can  be  estimated  with  use  of  a boundary  layer/loss 
calculation.  However,  to  rate  a set  of  similar  cascades,  a 
reasonable  criterion  can  be  constructed  with  the  maxima  and 
mimima  of  a blade  static  pressure  distribution  (Reference  4) . 

Based  on  the  parameter  of  Reference  4,  Equation  (2)  yields 
a value  to  compai'e  the  performance  of  a cascade  with  similar 
cascades. 


^®total  ~ I principal 
blade 


+ As 


splitter 

vane 


(2) 


where  As 


(P 

max 


P . ) 
min 


suction 

surface 


+ 


^1 


(P 

max 


P . ) 
min 


pressure 

surface 


In  order  to  determine  if  this  criterion  is  reasonable,  the  data 
presented  in  References  3 and  7 is  utilized  to  plot  the  experi- 
mental loss  coefficient,  versus  As.  This  is  shown  in  Figure 
7.  As  can  be  observed,  the  data  is  well  correlated  with  an 
increasing  As  yielding  an  increasing  loss.  Even  the  intro- 
duction of  a splitter  does  not  seriously  alter  the  correlation. 
Thus,  this  parameter  can  be  utilized  to  correlate  cascade  losses 
for  the  configurations  investigated. 


12 
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The  turning  performance  represents  a separate  problem.  An 
inviscid  calculation  is  not  complete  without  a specification  of 
total  turning.  This  is  often  done  by  imposing  a viscous  condi- 
tion on  the  trailing  edge  called  the  Kutta  condition  (see  Refer- 
ence 4) . Simply  stated,  the  Kutta  condition  prescribes  that  the 
stagnation  streamline  must  attach  to  a sharp  trailing  edge. 

This  condition,  along  with  the  requirement  of  periodicity,  is 
sufficient  to  establish  an  exit  flow  angle  for  a cascade  having 
subsonic  flow  at  the  exit.  The  practical  method  of  establishing 
the  validity  of  this  varies,  depending  on  what  is  viewed  to  be 
convenient  or  accurate  by  the  user  of  an  inviscid  program. 
Excellent  results  are  obtained  when  the  cascades  have  sharp 
trailing  edges  and  moderate  boundary  layers.  Unfortunately, 
the  cascade  of  this  study  violates,  to  some  degree,  both  of 
these  criteria.  The  result  is  that  the  deviation  angle  prediction 
should  be  less  accurate  than  that  noted  in  Reference  4. 

The  majority  of  the  calculations  for  the  transonic  relaxa- 
tion program  utilized  the  following  exit  boundary  conditions: 

o The  stagnation  point  on  the  trailing  edge  is 
specified. 

o The  shape  of  the  stagnation  streamline  is  adjusted 
periodically  throughout  the  calculation  procedure 
in  an  attempt  zo  enforce  periodicity. 

o A far-downstream  flow  angle  is  specified. 

Stagnation  streamlines  are  forced  to  this  angle 
at  a selected  distance  from  the  trailing  edge. 

In  light  of  the  uniqueness  of  turning  with  a specif i(?d.  stagnation 
point,  it  is  clear  that  the  boundary  conditions  are  over- 
specified.  Unless  the  specified  far-downstream  angle  is  the 
unique  angle  for  this  cascade,  periodicity  is  not  obtainable  no 
matter  how  much  adjustment  is  attempted.  Consequently,  some 
iteration  is  required  to  determine  the  proper  inviscid  turning 
for  a particular  cascade. 

If  it  is  desired  to  match  data  other  than  that  yielding  an 
ideal  Kutta  condition,  the  stagnation  point  can  be  moved  from  the 
trailing  edge  to  either  the  suction  or  pressure  surfaces.  How- 
ever, excess  movement  of  this  type  will  cause  an  irregularity  in 
the  grid  system  and  hence,  numerical  difficulty.  Consequently, 
turning  has  been  restricted  to  ideal  Kutta  conditions. 
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2.  PHASE  II  - ANALYSIS  OF  CASCADES  WITH  SPLITTER  VANES 

a.  Objectives 

The  objectives  of  this  phase  of  the  program  were  to  modify 
the  transonic-relaxation  program  to  accept  a splitter  of  arbi- 
trary geometry  and  to  compare  the  splitter  analysis  with  experi- 
mental data.  The  experimental  data  employed  in  this  comparison 
was  obtained  from  Reference  7.  In  the  test  program,  the  setting 
angle  of  the  splitter  vanes  could  be  altered  plus  or  minus  3 
degrees.  Therefore,  the  analysis  was  expanded  to  include  an 
investigation  of  the  effect  of  splitter  rotation  on  calculated 
splitter  loadings. 

b.  Cascade  Geometry  and  Test  Data 

Geometrical  details  of  the  main  blades  are  summarized  in 
Section  II.l.c.  Splitter  coordinates  chosen  to  numerically 
represent  the  blade  element  are  shown  in  Figure  8.  This  figure 
also  contains  pertinent  physical  characteristics  for  the 
splitter  vane. 

Experimental  data  for  the  ARL  supersonic  cascade  with 
splitters  is  presented  in  Figures  9,  10,  and  11  in  terms  of  wall 
static  pressure  versus  axial  distance.  The  cases  addressed  in 
these  figures  are: 

o A nominal  vane  setting  angle  of  39.75® 

o Rotated  by  plus  3®  or  a vane  setting  angle  of  42.75® 

o Rotated  by  minus  3®  or  a vane  setting  angle  or  36.75® 

When  the  cascade  was  operated  with  splitters,  the  vane  was 
instrumented  with  pressure  taps  to  obtain  a limited  number  of 
surface  static-pressure  measurements.  Information  such  as  main 
blade  and  splitter  static-pressure  is  included  in  these  figures. 

c.  Program  Modifications  for  Splitter  Vanes 

The  basic  program  for  computing  inviscid  transonic  flows  in 
a cascade  is  described  in  Section  Il.l.b.  That  program  was 
modified  to  incorporate  splitter  vane  geometry.  A total  of  22 
coordinate  points  on  each  surface  were  selected  to  represent  the 
splitter.  These  are  input  to  the  program  as  surface  boundary 
points . 

The  splitter  vane  is  handled  like  any  other  solid  wall  in 
that  the  normal  component  of  fluid  velocity  is  zero,  or  the 
derivative  of  the  potential  in  the  normal  direction  is  zero.  The 
coordinate  system  downstream  of  the  splitter  is  set  up  as  a 
dividing  streamline  with  no  flow  across  the  surface.  Initially, 
this  streamline  is  assumed  to  have  a linear  distribution  of  flow 
angle  from  the  trailing  edge  to  the  far  downstream  angle.  With 
each  major  iteration,  the  curvature  along  the  dividing  stream- 
line is  adjusted  to  equalize  the  fluid  velocities  on  each  side. 
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Figure  9.  Experimental  Data  for  Cascade  With 
Splitter  Vanes,  0®  Rotation. 
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In  essence,  the  two  passages  above  and  below  the  splitter  Vcuie 
are  isolated  and  treated  separately  beyond  the  splitter  leading- 
edge  station  except  for  this  curvature  correction  to  the  dividing 
streamline.  The  exit  grid  pattern  is  often  uneven  because  the 
total  number  of  grid  rows  is  constant  and  grid  spacing  changes 
between  passages.  However,  this  does  not  pose  a problem,  since 
the  far  downstrecun  exit  conditions  include  essentially  co  stant 
flow  angles  and  velocity  at  each  grid  intersection. 

A typical  grid  pattern  for  the  transonic  cascade  with 
splitters  is  shown  in  Figure  12,  Difficulty  was  experienced  in 
defining  the  grid  in  the  vicinity  of  the  splitter  leading  edge 
due  to  the  interpolation  scheme  for  locating  the  first  upper 
boundary  point  on  the  splitter  surface.  Under  certain  circum- 
stances, the  intersection  occurs  upstream  of  the  leading  edge 
station  and  the  prograun  fails.  This  is  remedied  by  changing  the 
upstreaun  grid  spacing  to  remove  the  upper  grid  intersection  from 
the  vicinity  of  the  splitter  leading  edge. 

Storage  schemes  had  to  be  modified  throughout  the  program 
because  the  location  of  the  leading  edge  in  the  grid  system  was 
unknown.  In  the  conventional  case,  odd  and  even  rows  of  grid 
points  are  different.  Even  rows  do  not  have  surface  intercepts, 
and  have  one  less  point  than  odd  rows,  as  illustrated  in 
Figure  13a.  However,  when  a splitter  is  introduced,  surface 
intercepts  may  occur  on  either  odd  rows  (Figure  13b)  or  even 
rows  (Figure  13c).  As  the  calculation  proceeds,  the  grid  system 
is  periodically  updated.  At  one  of  these  updates,  splitter 
intercept  patterns  may  change  from  even  to  odd  rows.  This 
causes  an  unavoidable  discontinuity  in  what  would  otherwise  be 
a suitable  iteration  process.  If  splitter  intercepts  are  on  an 
even  row,  exit  boundary  conditions  must  also  be  modified  since 
the  last  row  is  odd  but  does  not  contain  a splitter  surface 
intercept  as  is  the  case  with  standard  grid  systems. 

d.  Area  Calculations 

Area  calculations  were  performed  to  aid  in  uhe  interpreta- 
tion of  the  results  obtained  with  the  inviscid  transonic  pzogreun. 
Areas  were  obtained  geometrically  by  finding  a circle  tangpnt  to 
the  blade  surfaces.  The  location  of  the  center  of  this  circle 
was  used  to  obtain  a b-width.  The  area  was  then  determined  by 
applying  Equation  (3). 


A - 2rb  (3) 

This  process  involves  some  inaccuracy.  Fir-t,  the  flow  direction 
can  be  approximated  by  assuming  that  one-half  of  it  is  parallel 
to  one  surface  and,  the  other  half  is  parallel  to  ths  remaining 
surface.  The  process  also  assumes  that  a proper  average  b-width 
can  be  assigned  by  utilizing  a single  value  correlated  to  the 
circle  center.  The  level  of  this  inaccuracy  is  clearly  visible 
in  the  results  presented  in  Section  III. 
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Figure  13b.  Intercept  of 
on  Odd  Row. 


399-11 


3.  PHASES  III  AND  IV  - SPLITTER  REDESIGN 

A useful  approach  for  constructing  new  candidate  geometries 
was  to  begin  with  a solution  for  the  cascade  without  a splitter. 
Since  choking  is  a major  problem,  all  candidate  configurations 
must  be  selected  so  as  to  provide  an  adequate  area  distribution. 
At  the  same  time,  incidence  at  the  splitter  leading  edge  can  be 
estimated  from  free  streamlines.  To  obtain  free  streamlines,  a 
program  was  written  to  utilize  the  data  stored  on  disk  file  that 
resulted  from  the  inviscid  transonic  calculation.  This  program 
integrated  flow  from  surface  to  surface  along  successive  rows  of 
data.  The  location  of  a streamline  was  then  determined  and 
velocities  were  interpolated  to  this  streamline.  Pressure  and 
suction  surface  shapns  are  represented  by  well  reasoned  devia- 
tions from  the  streamline  shape.  This  technique  was  applied 
along  with  the  numerical  fluid  dynamic  solutions  discussed  in 
Sections  I and  II  to  redesign  splitter  configurations. 
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SECTION  III 


RESULTS 


1.  PHASE  I - ANALYSIS  OF  CASCADES  WITHOUT  SPLITTERS 

a.  Calculation  Method  - Surface  Pressures 

The  analyses  performed  in  this  investigation  assume  an 
inviscid  flow  field.  In  order  to  compare  measurements  of  blade 
surface  pressures  with  the  calculations,  a Mach  number  is  calcu- 
lated for  the  blade  surface  pressure  from  the  following  isen- 
tropic  relation;  1 


Y-1 

-| 

_2_ 

.J 

Y-1 

" 1 

. 

_ 

where  P is  the  local  static  pressure  on  the  blade  surface  and 
is  the  upstream  total  pressure. 

The  utilization  of  this  ideal  Mach  number  rests  on  two 
assumptions. 

(1)  A core  flow  with  essentially  no  total  pressure 

loss 

(2)  A boundary  layer  with  no  appreciable  normal  pres- 
sure gradients 

As  regards  the  first  assuir.ption,  cascade  data  taken  for  this 
configuration  indicates  that  considerable  loss  in  core  flow 
total  pressure  occurs  at  most  static  pressure  ratios.  The 
second  assumption  is  probably  true  in  unseparated  regions.  No 
experimental  evidence  exists,  however,  to  positively  establish 
the  size  and  extent  of  separations  on  the  blade  surface.  It  can 
be  inferred  that  they  must  exist  by  examining  both  ideal  and 
experimental  levels  of  diffusion. 

To  approximately  correct  the  value  of  Mach  number  computed 
from  Equation  4 for  total  pressure  loss,  the  following  expression 
is  a] so  utilized. 


1/2 


(5) 
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where 

Pt  - (C^-X)  CPji  - Pt2>^=X  '“t2 

and  where 

and  P^2  the  experimental  inlet 

and  exit  total  pressures  of  the  core  flow 


Several  subsequent  figures  show  both  methods  of  converting 
experimental  pressure  distributions  into  equivalent  Mach  numbe  ;s 
for  comparison  to  inviscid  calculations. 

b.  Comparison  of  Analysis  and  Test  Data 

It  is  difficult  to  compare  inviscid  flow  calculations  with 
experimental  results  w?.ere  viscous  effects  are  so  evident.  It 
is  not  possible  to  make  a single  comparison  of  all  the  important 
parameters  across  the  transonic  cascade.  The  following  table 
summarizes  the  problem  and  indicates  how  the  final  comparisons 
are  made. 


TABLE  II.  PARAMETiPS  SELECTED  FOR  PHASE  I ANALYSIS. 


Analytical  Data 

Experimental  Data 

1 

0 

Turn- 

ing 

^2 

AVR-Axial 

Velocity- 

Density 

Ratio 

Static 

Pres- 

sure 

Ratio 

0 

Turn- 

ing 

®2 

AVR-Axial 

Velocity- 

Density 

Ratio 

Static 

Pres- 

sure 

Ratio 

34.65 

32.20 

2.08 

2.72 

— 

— 

— 

— 

m 

42.35 

24.00 

2.53 

2.36 

36.90 

29.95 

2.07 

2.038 

0.808 

45.05 

21.80 

2.62 

2.26 

37.39 

79.47 

2.15 

1.968 

0.9 

46.30 

20.55 

2.67 

2.05 

37.25 

29.60 

2.24 

1.821 

Initial  conditions  for  the  analysis  were  fixed  at  an  inlet 
Mach  number  of  1.46  and  an  inlf t flow  angle  of  66.85  degrees  as 
in  the  experimental  setup.  Analyses  were  than  conducted  for 
three  exit  Mach  numbers;  the  latter  two  cases  (M2  = 0.76  and  0.9) 

corresponded  to  experimental  conditions.  The  0.6  Mach  number 
case  corresponds  to  inviscid  flow  with  no  aerodyneimic  blockage. 
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This  case  has  an  axial  velocity-density  ratio  similar  to  the 
experimental  condition  at  M = 0.76  and,  therefore,  provides  an 
additional  means  for  comparing  the  two  results.  As  indicated  in 
Table  II,  the  turning,  axial  velocity-density  ratio,  and  static 
pressure  ratio  are  different  when  analytical  and  experimental 
results  are  compared  at  the  same  exit  Mach  number.  However,  by 
comparing  cases  with  similar  axial  velocity-density  ratios,  the 
turning  is  more  nearly  the  same  but  exit  Mach  numbers  and  static 
pressure  ratios  are  different.  Thus,  the  former  comparison 
Cequal  M2)  is  employed  to  match  analytical  and  experimental  blade 
loadings  while  the  latter  comparison  (Equal  AVR)  is  made  to 
evaluate  the  flow  turning  conditions.  This  is  illustrated  in 
Figure  14. 

Analytical  results  for  exit  Mach  numbers  of  0.6,  0.76, 

0.808,  and  0.9  are  shown  on  the  subparts  of  Figures  15  through 
18  respectively.  These  figures  are  divided  into:  (a)  surface 

Mach  number,  (b)  the  grid  system,  (c)  the  variation  of  f4ach 
numbers  along  quasi-streamlines  and,  except  for  Figure  15, 

(d)  contours  of  constant  Mach  number.  Surface  Mach  numbers  are 
calculated  directly  with  use  of  the  program  described  in  Section 
II.  Experimental  values  are  calculated  using  one  of  the  two 
methods  described  in  (a)  above.  The  grid  system  shows  Mach 
lines  in  the  supersonic  region.  This  permits  interpretation  of 
the  effects  of  reflections  of  one  surface  on  another.  The 
streamline  curves  represent  an  exploded  plot  of  Mach  number  along 
lines  obtained  by  connecting  grid  points  in  a streamwise  fashion. 
The  number  imprinted  at  the  start  of  each  line  is  an  adder  used 
to  spread  the  curves.  The  Mach  number  read  is  equal  to  the  true 
Mach  number  plus  the  adder.  The  last  plot  is  a plot  of  lines  of 
constant  Mach  number  within  the  cascade.  No  experimental  Mach 
numbers  are  shown  for  the  0.6  Mach  number  case  since  it  repre- 
sents the  b-width  distribution  of  the  cascade  end  wall  without 
any  aerodynamic  blockage  effects. 

A comparison  with  time-dependent  methods  was  made  to  verify 
the  accuracy  of  the  relaxation  method.  The  surface  Mach  number 
distributions  resulting  from  this  comparison  are  shown  in 
Figure  19.  Two  different  b-width  distributions  are  shown.  As 
one  of  its  boundary  conditions,  the  time-dependent  calculation 
always  satisfies  an  ideal  Kutta  condition.  The  exit  Mach  number, 
however,  will  vary  with  the  downstream  flow  angle  required  to 
meet  the  Kutta  condition.  Thus,  to  precisely  match  data  when 
using  the  time-dependent  method,  a variety  of  axial  velocity- 
density  ratios  must  be  applied  until  one  is  found  that  yields  the 
exit  Mach  number  desired.  Table  III  summarizes  the  results  of 
the  comparison  of  the  two  calculation  methods  with  test  data. 

The  numerical  smoothing  of  the  time-dependent  method  provides  a 
much  smoother  surface  Mach  number  distribution  than  that  of  the 
relaxation  solution.  However,  the  general  features  of  the  relax- 
ation process  are  all  confirmed.  In  particular: 


AXIAL  VELOCITY  - DENSITY  RATIO 


0.5  0.6  0.7  0.8  0.9  1.0  1.1 


M2  ~ EXIT  MACH  NUMBER 


Figure  14.  Inviscid  Simulation  of  Real  Fluid 
Conditions  in  a Cascade. 
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Figure  15a.  Surface  Mach  Number  Distribution  for 
Compressor  Cascade  at  M-  of  0.6. 
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Figure  16a.  Con^>arison  of  Experimental  Blade 
Surface  Mach  Number  Distribution 
with  Results  of  Calculation  for 
M2  “ 0.760. 
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Figure  17a. 


Comparison  of  Experimental  Blade 
Surface  Mach  Number  Distribution 
with  Results  of  Calculation  at 
M2  of  0.808. 
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Figure  18  a.  Comparison  of  Experimental  Blade  Surface 


Mach  Number  Distribution  with  Results  of 
Calculation  at  M2  of  0.900. 


Figure  18b.  Non-Orthogona 
Cascade  at  M- 
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(1)  The  Mach  numbers  at  the  leading  edge  pressure 
surface  are  similar  to  the  relaxation  results,  and  definitely 
higher  than  experimental  values. 

(2)  The  Mach  numbers  near  the  trailing  edge  pressure 
surface  are  significantly  below  experimental  values,  but  agree 
with  relaxation  calculations. 

(3)  Suction  surface  Mach  number  distributions  over  the 
first  60  percent  of  the  blade  agree  reasonably  well  with  both 
relaxation  and  experimental  values. 

(4)  The  relaxation  solutions  indicate  a stronger  com- 
pression on  the  latter  40  percent  of  the  suction  surface  than 
does  the  time-dependent  method.  However,  the  time-dependent 
method  is  still  stronger  than  the  experimental  values. 

In  addition,  time-dependent  methods  predict  larger  turning  than 
actually  occurs  although  it  is  not  quite  as  large  as  that  indi- 
cated by  relaxation  solutions.  See  Table  III. 

c.  Effect  of  b-Width  Variation 


In  order  for  the  computations  to  achieve  the  measured  outlet 
Mach  number  for  specified  values  of  inlet  Mach  number  and  flow 
angle,  the  effective  outlet  b-width  must  be  less  than  the  geome- 
trical b-width.  This  is  physically  a consequence  of  boundary 
layer  growth  through  the  cascade. 

Equation  (1)  is  the  continuity  equation  for  the  overall 
flow  between  the  inlet  Station  1 and  the  downstreeun  Station  2 . 

For  the  ratio  of  specific  heats,  Yf  equal  to  1.4  and  the  speci- 
fied inlet  flow  test  conditions  (M^^  = 1.460,  = 66.850)*, 

Equation  (1)  becomes 


0.197810 


(l  + 0.2M2)^ 

M2  cos  6^ 


(7) 


This  equation  shows  that  for  values  of  M2  specified  by  test  condi- 
tions, the  overall  effective  b-width  ratio  depends  upon  the  value 
of  $2  required  to  satisfy  the  Kutta  condition. 


*The  lack  of  periodicity  in  the  analyses  showed  that  periodic 
upstream  flow  was  not  obtained  with  these  conditions. 
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The  overall  change  in  the  effective  h-width  does  not  yield 
sufficient  information  to  determine  the  variation  of  the  b-width 
within  the  cascade.  Instead,  assumptions  are  required  concerning 
the  variation  of  the  b-width.  It  was  assumed  that  the  b-width 
variation  was  primarily  controlled  by  the  variation  of  the  geome- 
trical b-width  through  the  cascade.  Then,  the  change  in  local 
b-width  from  its  geometrical  value  can  be  expressed  by 


Ab(X) 


where  f(X)  is  some  function  of  the  axial  coordinate,  X,  through 
the  cascade  expressing  the  local  change  in  b-width  in  relation  to 
the  overall  change  in  b-width. 

The  following  three  assumptions  for  f(X)  were  investigated: 

Assumption  1:  f(X)  = a^  + a^^X 

with  boundary  conditions 

o f(X)  = 0 at  X = -0.0546  (0.1-inch  axially  upstream 

of  leading  edge) 

o f(X)  = 1 at  X = 1.1458  (approximately  0.5%  of 

forward  of  the  region  of  contraction 
of  cascade  end  walls) 

Assumption  2:  f (X)  = 

with  boundary  conditions 

o f(X)  = 0 at  X = 0.5 

o f(X)  = 1 at  X = 1.1458 

2 3 

Assumption  3t  f(X)  “ + a2X  + a^X 

with  bouncary  conditions 

o f(X)  = 0 at  X = -0.0546 
o df(X)/dX  = 0 at  X = -0.0546 
o f (X)  = 1 at  X ='  1.1458  Cy. 
o df(X)/dX  = 0 at  X = 1.1458  C„ 

The  variation  of  f(X)  is  illustrated  in  Figure  20  for  these 
three  assumptions. 
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Figure  20.  Curves  Describing  f(X). 
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Excunples  of  the  blade  surface  Mach  number  distributions 
resulting  from  the  three  assumptions  of  the  functional  form  of 
f(X)  are  illustrated  in  Figure  21.  The  results  are  for  the  out- 
let flow  angles  that  approximately  satisfy  the  Kutta  condition. 
Significant  differences  are  noted  in  the  trailing  edge  region  and 
on  the  pressure  surface  from  X/C^  = 0.5.  The  calculation  results 
yield  different  values  for  the  frow  angle  at  the  trailing  edge  as 
indicated  in  Table  IV.  As  may  be  seen,  the  deviation  angle  at 
the  trailing  edge  varies  from  near  zero  to  2.5  degrees,  depending 
upon  the  manner  in  which  f(X)  is  represented. 

TABLE  IV.  ESTIMATES  OF  FLOW  ANGLE  AT  THE 
TRAILING  EDGE  PLANE  UNDER  THE 
INFLUENCE  OF  b-WIDTH  DISTRIBUTION. 


Description 

^2 

«2 

®TE 

Experimental  results 
(P2/P]^  = 1.821) 

29.600® 

0.900 

- 

22.534®* 

Linear  f(X)  from  X/C„  = 
-0.0546  (Assumption'll) 

20. 55® 

0.9514 

0.8766 

22.  3® 

Linear  f(X)  from  X/C^  = 
0.5  (Assumption  2) 

21.0® 

0.9581 

0.8281 

24.5® 

Smooth  f(  . (Assumption  3) 

25.0® 

0.9617 

0.9597 

25.1® 

*'i railing  edge  blade  angle  on  mean  camber  line  (see  Table  I)  . 


The  choice  between  the  three  assumed  functions  f(X),  as  well 
as  any  other  '•''.at  may  be  chosen,  is  purely  arbitrary  since  no 
means  exist  r determining  the  variation  of  the  effective 
b“Width.  Fui  .:hermore,  comparison  of  the  calculated  blade  surface 
Mach  number  distributions  with  measured  values  is  not  xikely  to 
confirm  a choice  since  the  measured  values  necessarily  include 
the  viscous  effects  of  the  real  flow  (i.e.,  blade  boundary  layer 
growth  due  to  shock-boundary  layer  interaction — an  effect  that 
leads  to  smoothing  the  blade  surface  pressure  distributions  in 
comparison  to  inviscid  flow  calculations) . For  this  reason,  the 
simplest  assumption  (Assumption  1)  was  'ised  in  this  investigation. 
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Figure  21.  Comparison  of  Blade  Surface  Mach  Number  Distributions 
for  Various  Assumptions  on  the  Behavior  of  f(X). 
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d.  Pressure  Surface  Leading  Edge 

On  both  a qualitative  and  quantative  basis,  the  Mach  number 
distribution  near  the  pressure  surface  leading  edge  differs 
markedly  from  that  of  the  data.  All  relaxation  calculations 
show  a compression  over  the  first  10  percent  of  the  blade 
followed  by  a sharp  but  brief  re-acceleretion  and  a strong  shock 
wave.  Test  data  shows  a compression  near  the  leading  edge  that 
is  a strong  function  of  cascade  back  pressure  without  any  re- 
acceleration  such  as  that  observed  by  the  calculation.  The  time-  I 

dependent  calculation  falls  between  the  data  and  the  relaxation  | 

calculation,  indicating  a weaker  compression  b”t  no  re- 
acceleration. 

Analysis  of  the  area  distribution  through  the  cascade  indi- 
cates that  the  area  continuously  decreases  from  the  leading  edge 
until  a minimum  is  reached  midway  through  the  Cc.scade.  This 
minimum  has  approximately  8 percent  less  area  than  exists  at  , 

the  leading  edge.  Figure  22  shows  the  approximate  location  of 
this  minimum.  The  minimum  lies  between  the  pressure  surface  at  | 

approximately  0.31  axial  chord  and  the  suction  surface  at  0.72  ■ 

axial  chord.  The  calculated  strong  compression  occurs  just  i 

downstream  of  this  location.  Correspondingly,  the  strong  com- 
pression evidenced  on  the  suction  surface  occurs  very  near  0.72  i 

of  axial  chord.  The  presence  of  a strong,  nearly  normal  wave  in  | 

this  region  (shown  in  Figures  16d,  17d,  and  18d)  seems  to  be  1 

related  to  this  minimum  passage  area.  The  lack  of  a similar  j 

strong  compression  in  the  data  is  evidence  of  strong  viscous  I 

action  in  this  region.  \ 

A comparison  of  Mach  numbers  calculated  by  relaxation  and  ! 

derived  from  experimental  data  is  chown  in  Figure  23.  The  I 

oblique  shocks  going  from  upstream  conditions  (M  = 1.460,  s 

8 = 66.85”)  to  the  local  surface  angle  are  shown  in  Figure  23  to  | 

permit  x reasonable  compression  value  to  be  estimated.  Two  i 

branches,  a strong  obJique  and  weak  oblique,  are  indicated.  The  | 

relaxation  analysis  shows  compressior « -.'sbly  consistent  ] 

with  the  weak  oblique  shock  relations  -.c  ..  7 maximum  static 
pressure  ratio,  the  data  yields  a substan.  . \ stronger  com- 
pression than  either  the  relaxation  analys;'.s  or  the  oblique  shock 
relations.  A compression  of  this  magnitude  wouxd  require  an 
oblique  shock  yielding  in  excess  of  10.5  degrees  of  deflection. 

This  is  5 degrees  of  deflection  past  the  local  metal  angle.  Such 
a deflection  is  possible,  since  the  static  pressure  ratio  exceeds 
that  needed  to  separate  a boundary  layer  (static  pressure  from 
far  upstream  to  the  first  point  measured  on  the  pressure  surface 
for  the  M2  = 0.76  case  is  1.78).  However,  other  effects  could 
also  be  responsible  since  this  case  is  far  more  complex  than 
those  usually  treated  by  simple  oblique  shock  relations.  Waves 
generated  abreast  of  the  blade  should  be  weak  since  the  unique 
inlet  flow  angle  is  nearly  the  suction  surface  angle.  The 
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relatively  strong  contraction  in  b-width  complicates  the  under- 
standing of  the  phenomenon  in  the  leading  edge  region.  For 
example,  an  unloaded  streamline  should  have  nearly  constant  Mach 
number  along  it.  When  there  is  no  change  of  b-’width,  such  a 
streamline  is  straight.  When  b-width  is  not  constant,  the  follow- 
ing relations  can  be  applied: 

p V cosB  b = p,V  cosB  b^  Continuity  (5) 

^ ^ of  flow 


V sinB  = sing.  Constant 

tangential  (7) 

momentum 


The  result  is  that  an  unloaded  streamline  will  have  a flow  angle 
given  by  Equation  (8) 


tanS  = 


Pb 

Plbi 


tang^ 


(8) 


Often,  when  8 is  large,  the  density  ratio  is  near  unity  and 
Equation  (7)  becomes  dominantly  geometric.  Thus,  the  suction 
surface  even  though  curved,  provides  almost  constant  velocity 
near  the  leading  edge.  If  the  pressure  surface  has  a radius  of 
curvature  less  than  that  of  a free  streamline,  it  will  be  com- 
pressive (curve  b of  Figure  24) . If  it  has  a radius  of  curvature 
greater  than  a free  streamline,  it  will  be  expansive  (curve  C of 
Figure  24) . A streamline  approaching  from  far  upstream  must  he 
compressive  as  illustrated  by  Figure  25  (note  the  suction  surface 
lies  fairly  close  to  a free  streamline  and  is  only  slightly 
expansive) . If  the  density  variations  are  neglected.  Equation 
(8)  can  be  differentiated  simply  to  yield  Equation  (9) . 

d(tan8)  _ ^^"^1  d(b) 

d(X/Cjj)  - -b]“  r(Vcp"  (9) 


This  equation  can  be  used  to  estimate  whether  a surface  should 
be  compressive  or  expansive.  A plot  of  Equation  C9)  compared  to 
the  actual  derivative  of  tan8  on  the  pressure  surface  is  shown  in 
Figure  26.  Note  that  if  the  flow  upstream  were  aligned  with  the 
pressure  surface,  it  would  expand  along  the  pressure  surface  until 
a point  in  excess  of  20  percent  of  chord  is  reached,  whereupon  it 
would  compress.  This  indicates  that  it  is  reasonable  to  expect 
the  following  on  the  pressure  surface. 


COMPRESSIVE  (a)  expansive 

STREAMLINE  FREE  STREAMLINE 

STREAMLINE 
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Figure  25.  Flow  Angles  Along  Suction  Surface, 
Pressure  Surface,  and  a Free 
Streamline, 
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Figure  26.  Derivative  of  Pressure  Surface  Angle  Versus 
Axia;  Position. 
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(1)  Compression  from  the  upstream  direction  to  the 
local  flow  direction  at  leading  edge 

(2)  Expansion  behind  the  compression  due  to  the 
sharp  b-width  reduction  on  a surface  of  low  curvature 

(3)  Subsequent  compression  due  both  to  increased 
surface  curvature  and  a decrease  in  available  flow  area 

This  general  structure  is  exactly  what  was  computed.  Its  absence 
in  the  data  is  attributed  to  strong  viscous  effects  in  the  lead- 
ing edge  region.  The  effect  on  design  could  be  significant 
because  the  inviscid  calculation  predicts  high  Mach  numbers  imme- 
diately ahead  of  the  splitter.  The  data  indicates  a much  lower 
level  of  Mach  number,  at  least  on  the  pressure  surface. 

e.  Exit  Region 

It  has  been  established  experimentally  that  an  important 
parameter  for  determining  exit  deviation  angle  is  axial  velocity 
density  ratio.  At  the  same  time,  exit  deviation  angle  is  largely 
independent  of  Mach  number  when  large  viscous  effects  are  not 
present.  Consequently,  to  predict  turning,  axial  velocity  den- 
sity ratio  should  be  matched.  When  losses  are  high,  an  inviscid 
calculation  cannot  match  the  axial  velocity-density  ratio  and  the 
static  pressure  rise  of  the  cascade  even  approximately.  Adjust- 
ment of  b-width  will  reproduce  either  exit  Mach  number  or  axial 
velocity-density  ratio,  but  not  both.  Consequently,  when  b-width 
distributions  are  adjusted  to  match  exit  Mach  number,  the  calcu- 
lated turning  is  greater  than  is  observed  experimentally.  The 
result  is  more  turning  near  the  trailing  edge  and  rather  large 
discrepancies  in  loading  distribution  in  the  latter  portion  of 
the  blade. 

f.  Choking 

It  is  important,  when  considering  the  application  of  split- 
ters, to  determine  the  ability  of  the  prediction  method  to  pre- 
dict choke.  Boundary  conditions  for  the  inviscid  transonic 
releixation  are  such  that  a constant  flow  rate  is  forced  into  the 
xcascade.  If  this  exceeds  the  flow  capacity,  a "spill"  in  the 
choked  region  will  occur  which  causes  the  coefficient  matrix  to 
become  singular  and  the  relaxation  solution  to  diverge.  Study  of 
several  cascades  indicates  that  this  "choking"  failure  takes 
place  when  the  flow  is  just  less  than  that  yielding  a local  A/A* 
of  unity. 
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For  the  cascade  under  study  in  Phase  I,  the  area  decreases 
to  a minimum  (approximately  8 percent  less  than  the  inlet) . 

Since  the  A/A*  for  an  inlet  Mach  number  of  1.460  is  1.1501  this 
is  not  sufficient  to  cause  choking.  The  area  then  increases  to 
the  trailing  edge.  The  minimum  passage  b-width  occurs  down- 
stream of  the  trailing  edge.  Figure  27  shows  the  combination  of 
turning  and  axial-velocity-density  ratio  that  results  in  choking 
(A/A*  =1.0)  at  this  downstream  station.  Relaxation  calculations 
were  possible  at  axial  velocity-density  ratios  up  to  at  least 
2.67.  Calculations  were  not  possible  when  axial  velocity-density 
ratios  were  increased  beyond  a certain  point.  During  the  calcu- 
lations to  determine  the  Kutta  condition,  results  were  obtained 
to  an  exit  Mach  number  of  0.98.  Normal  calculations  to  deter- 
mine the  Kutta  condition  consisted  of  varying  the  exit  angle 
with  a constant  axial  velocity  density  ratio.  Variations  below 
the  Kutta  line  (Figure  27)  were  observed  to  be  severely  restric- 
ted by  "choking"  failures. 

g.  Conclusions  - Phase  I 

A review  of  the  analysis  completed  in  Phase  I suggests  the 
following; 


(1)  The  transonic  relaxation  solution  represents  a 
reasonable  approximation  to  the  inviscid  flow  within  the  cascade. 

(2)  Substantial  deviations  between  experimental  data 
and  inviscid  calculations  were  observed.  Principal  deviations, 
from  a design  point  of  view,  are: 

o The  lack  of  a strong  pressure  surface  re- 
acceleration followed  by  compression.  The  invis- 
cid solution  implies  relatively  high  Mach  number 
abreast  of  the  location  of  a splitter  leading  edge 

o Larger  turnings  than  those  observed  at  a given 
exit  Mach  number 

o Less  diffusion  on  the  pressure  surface 

o The  exit  flow  angle  at  the  Kutta  condition  was 
appreciably  influenced  by  b-width  shape  and 
distribution 

(3)  There  is  general  agreement  between  the  experi- 
mental data  and  the  calculations  along  the  suction  surface. 

It  should  be  possible,  with  use  of  the  inviscid  solution,  to 
determine  the  configuration  offering  the  lowest  overall  surface 
diffusion.  Consequently,  the  lowest  loss  configuration  could 
also  be  determined.  However,  a reasonable  design  estimate  of 
deviation  angle  probably  cannot  be  based  on  the  inviscid  solution. 
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2.  PHASE  II  - ANALYSIS  OF  CASCADES  WITH  SPLITTER  VANES 

a.  Relaxation  Solution 

The  results  of  relaxation  calculations  with  the  ARL  splitter 
at  two  different  splitter  settings  are  presented  in  Figures  28a 
through  28d,  and  29a  through  29d.  The  following  information  is 
shown  in  each  figure. 

28a  and  29a  - Final  grid  system 

28b  and  29b  - Mach  number  distribution  along  quasi- 
streamlines 

28c  and  29c  - Lines  of  constant  Mach  number 

28d  and  29d  - Surface  static  pressures 

Relaxation  calculations  were  performed  with  a b-width  distribu- 
tion that  corresponds  to  the  cascade  test  metal  dimensions. 

Mach  number  distributions  are  for  the  highest  experimental 
static-pressure  ratio,  even  though  the  experimental  exit  Mach 
number  is  always  higher  than  the  0.6  exit  Mach  number  computed. 

All  calculations  with  the  exit  Mach  number  near  the  experi- 
mental values  indicated  that  the  cascade  was  choked.  These  cal- 
culations wrjre  performed  with  a natural  flow  split.  That  is,  the 
one  with  the  continuous  potential  across  the  downstream  boundary 
between  the  upper  and  lower  passages.  Attempts  were  made  to 
modify  the  calculated  flow  split  between  upper  and  lower  passages. 
This  could  be  done  by  increasing  the  potential  on  the  far  down- 
stream boundary  of  one  passage  relative  to  the  other.  However, 
it  would  also  result  in  higher  velocities  in  one  passage  than  it 
would  in  the  other.  This  complicates  the  calculation  of  the 
shape  of  the  boundary  between  passages  since  the  curvature  of 
this  slip  line  is  adjusted  to  yield  a static  pressure  baJance  on 
the  common  contour  between  passages  downstream  of  the  splitter. 

The  changes  required  to  acconplish  this  modified  flow  split 
correctly  were  extensive  and  outside  the  scope  of  the  present 
program.  The  observed  choking  led  directly  to  the  examination 
of  flow  area,  which  is  discussed  in  Paragraph  b. 

b.  Area  Calculation 

The  overall  area  distribution  through  the  passage,  with  and 
without  splitters  is  shown  in  Figure  30.  The  introduction  of  the 
splitter  results  in  a 2-percent  reduction  of  the  minimum  flow 
area.  However,  the  overall  area  is  well  above  the  15-percent  con- 
traction required  to  choke  the  cascade.  Individual  passage  areas 
yield  a quite  different  result.  An  estimate  of  the  flow  split 
can  be  obtained  by  observing  the  streamline  (without  the  splitter) 
passing  through  the  location  of  the  splitter  leading  edge.  This 
is  done  using  the  flow  solution  without  splitter  discussed  in 
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Section  II.  The  result,  given  in  Figure  31,  shows  approximately 
60  percent  of  the  flow  passing  between  the  splitter  leading  edge 
and  the  main  blade  suction  surface.  The  area  calculated  for  each 
streamtube  relative  to  A*  (the  choking  area)  of  the  upstream  flow 
is  shown  in  Figure  32.  The  net  error  in  the  geometric  approach 
to  area  calculation  is  evident,  since  A/A*  should  never  be  less 
than  unity.  This  is  true  since  the  area  comes  from  a solution 
that  actually  passed  the  flow.  Nonetheless,  the  result  is  clear. 
The  lower  passage  (between  the  splitter  pressure  surface  and  the 
main  blade  suction  surface)  is  extremely  deficient  in  area,  while 
very  large  areas  are  available  in  the  upper  passage.  An  analysis 
of  this  indicates  that  the  lower  passage  utilizes  a b-width  based 
on  an  average  X that  is  larger  than  the  upper  passage.  Since  b- 
width  converges  rapidly  with  X,  the  result  is  a choked-lower 
passage. 


c.  Discussion  of  the  Results 

The  effect  of  choking  the  lower  passage  is  markedly  differ- 
ent between  the  experimental  data  and  the  calculation.  In  the 
relaxation  calculation,  a solution  is  not  possible  once  the 
passage  becomes  choked.  In  the  experiment,  losses  rise  sharply 
in  the  choked  (lower)  passage  and  a flow  shift  towards  higher 
upper  passage  flow  results.  This  yields  an  incidence  on  the 
splitter  that  is  quite  different.  Except  at  the  highest  stagger, 
the  data  indicates  a positive  incidence  on  the  splitter.  Calcula- 
tions show  a slight  negative  incidence  at  design  stagger.  The 
incidence  shifts  to  positive  as  the  stagger  is  decreased,  indicat- 
ing a neutral  incidence  between  39.75-degrees  and  36.75-degrees 
splitter  stagger.  Comparison  of  the  data  indicates  that  experi- 
mental incidence  was  approximately  3-degrees  higher  than  calcu- 
lated incidence. 

Further  evidence  of  this  choking  can  be  obtained  from  the 
main  blade  pressure  distributions.  Pressure  surface  Mach  num- 
bers show  a sharp  rise  in  the  vicinity  of  the  splitter  leading 
edge,  which  is  not  present  when  the  splitter  is  removed.  This  is 
consistent  with  an  increase  in  flow  in  the  upper  passage  over 
what  would  normally  be  experienced.  Since  the  method  used  will 
not  calculate  any  case  beyond  the  choke  point,  none  of  these 
effects  are  visible  in  the  calculation.  An  attempt  was  made  to 
modify  the  flow  split  away  from  one  yielding  lines  of  constant 
potential  at  the  exit.  However,  the  coding  changes  required  to 
do  this  were  beyond  the  scope  of  the  program. 
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3.  PHASES  III  ANH  iV  - SPLITTER  REDESIGN 
a.  Methodology 

Based  on  the  results  of  Phases  I and  II,  the  following  con- 
clusions were  achieved. 

(1)  Viscous  effects  are  large  enough  to  reduce  the 
quality  of  agreement  between  the  calculation  and  the  experiment. 

(2)  The  current  cascade  has  a splitter  configuration 
that  (because  of  loss  and  boundary  layer  buildup)  appears  to 
operate  v/ith  a choked  lower  passage. 

(3)  The  choked  lower  passage  causes  a higher  than 
predicted  incidence  on  the  splitter  of  approximately  3 degrees. 

(4)  A comprehensive  comparison  of  experimental  data 
and  calculations  is  not  possible  with  tested  splitter  config- 
urations due  to  the  absence  of  calculations  beyond  choke. 

Even  though  comprehensive  comparisons  were  not  possible  in  Phase 
II,  a redesign  was  attempted  in  Phase  III.  The  choking  condition 
observed  in  Phase  II  must  be  detrimental  to  cascade  performance 
since  it  creates  high  losses  in  the  lower  passage  and  large  over- 
velocities on  the  splitter  suction  surface  due  to  the  induced 
positive  incidence.  If  a configuration  exists  that  possesses 
adequate  area  and  reasonable  splitter  pressure  distributions, 
its  performance  should  be  better. 

An  examination  of  solutions  without  splitters  (Figure  32) 
indicates  that  simply  moving  the  splitter  closer  or  further  away 
from  the  suction  surface  will  not  help.  The  alternatives  seem  to 
be  restricted  to  one  of  the  following: 

(1)  Increasing  splitter  length  until  larger  flow 
area  is  encountered 

(2)  Shortening  the  splitter  until  larger  flow  area 
is  encountered 

(3)  Changes  in  camber  line  shape 

(4)  Changing  b-width  distribution 

(5)  Changing  main  blade  shape 

Alternatives  (4)  and  (5)  are  changes  to  the  entire  rotor  and, 
thus,  are  outside  the  scope  of  feasible  changes.  Consequently, 
Alternatives  (1),  (2),  and  (3)  were  investigated.  Shortening  the 
splitter  is  of  more  academic  interest  than  practical  value,  since 
it  would  probably  yield  a rotor  with  stress  and  vibration 
problems.  Lengthening  the  splitter  would  move  the  splitter  lead- 
ing edge  into  high  Mach  number  regions,  which  would  increase  the 
loss  of  the  splitter  through  increased  dynamic  head. 
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Changes  in  camber  line  might  be  effective,  as  a result  of 
improved  passage  area,  except  that  the  splitter  is  very  close  to 
the  minimum  cascade  area. 

b.  Increased  Length 

It  would  be  desirable  to  place  the  splitter  on  a streamline 
at  a position  where  adequate  flow  area  would  occur  both  above  it 
and  below  it.  The  shape  of  the  pressure  surface  would  then  be 
adjusted  so  that  the  area  distributions  would  never  fall  below  a 
certain  minimum  value.  Such  a design  is  shown  in  Figure  33.  The 
leading  edge  is  located  at  the  55-percent  flow  streamline.  The 
pressure  surface  is  selected  to  provide  a constant  A/A*  of  1,05. 
The  resulting  area  distribution  is  shown  in  Figure  34.  The  grid 
system  is  shown  in  Figure  35.  Unfortunately,  this  configuration 
results  in  a large  negative  incidence.  The  calculation  indicates 
a large  expansion  in  the  lower  passage,  rising  Mach  number  near 
the  splitter  leading  edge,  and  a choked  lower  passage.  This 
choking  occurs  even  though  on  a one-dimensional  basis,  sufficient 
area  exists  to  pass  the  flow.  The  expansion  was  strong  enough 
to  dominate  the  flow  field  and  raised  the  Mach  number  to  a point 
where  the  lower  passage  effective  A/A*  was  less  than  the  actual 
geometric  A/A*  available.  This  could  result  in  a normal  shock 
in  front  of  the  splitter  along  with  a flow  shift  toward  the  upper 
passage  in  an  actual  cascade  test.  Clearly,  this  would  be  an 
undesirable  condition. 

To  avoid  this,  the  leading  edge  would  need  to  be  moved  for- 
ward and  the  camber  line  established  in  a more  gradual  S-shape  to 
control  the  expansion  necessary  for  a contour  with  sufficient 
area  to  avoid  choking.  This  problem  is  complicated  by  the 
decreasing  b-width  which  serves  to  increase  expansion  levels  fur- 
ther than  would  be  the  case  with  a two-dimensional  cascade.  As 
a result  of  this  study,  it  did  not  appear  that  increased  splitter 
length  would  yield  a practical  solution  even  with  a change  in 
splitter  shape. 

c.  Shortened  Splitter 

A somewhat  shortened  splitter  is  illustrated  in  Figure  36. 

An  area  distribution  is  shown  in  Figure  37.  Examination  of 
Figure  37  indicates  that  the  lower  passage  area  distribution  is 
relatively  flat  up  to  the  trailing  edge.  Consequently,  for  a 
large  improvement,  an  extremely  shortened  splitter  would  be 
needed.  Little  improvement  would  be  expected  with  the  shortened 
splitter  shown  in  Figure  36,  Figures  38a  through  38d  show  the 
results  of  a calculation  with  no  aerodynamic  blockage.  More 
severe  b-width  distributions  yielded  a choked  lower  passage  as 
was  the  case  with  the  design  splitter.  Further  efforts  along 
this  line  were  abandoned  since  the  splitter  would  need  to  be 
impracticably  short  to  avoid  a lower  passage  choke. 
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Figure  33.  Modified  Splitter  for  Cascade 
(Increased  Length) . 
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Figure  34.  Area  Distributions  for  Modified  Splitter 
Passages  (Increased  Length) . 
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Figure  38d.  Surface  Mach  Number  Distribut 
With  Modified  Splitter  (Short 
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SECTION  IV 
CONCLUSIONS 


An  inviacid  calculation  can  be  used  to  analyze  fluid  dynamic 
elements  when  viscous  effects  are  sufficiently  small  to  permit 
the  flow  to  be  dominated  by  pressure  forces.  When  visions 
effects  are  veary  low,  good  calculated  predictions  of  pressure 
distributions  can  be  obtained.  This  has  been  observed  to  be  the 
case  for  most  low-speed  cascade  data.  As  viscous  effects  become 
more  prominent,  the  inviscid  calculation  does  not  accurately  pre- 
dict the  effects  of  boundary  layer  growth  and  increased  losses  on 
the  pressure  distributions  within  the  cascade.  However,  even 
when  v:‘ scous  effects  are  large,  an  invsicid  calculation  can  be 
used  to  identify  regions  of  high  loading  and  to  provide  an  in- 
sight to  the  major  sources  of  high  losses. 

When  a splitter  vane  is  added  to  the  cascade,  the  system, 
becomes  more  complex.  Viscous  effects  that  caused  only  modest 
disagreements  between  calculated  and  experimental  velocity  dis- 
tributions in  the  plain  cascade  are  now  largely  responsible  for 
establishing  the  flow  capacity  of  the  passages  on  each  side  of 
the  splitter.  Areas  of  high  loss  reduce  the  local  flow  capacity 
and  cause  flow  shifts  relative  to  the  splitter  that  are  not  pre- 
dicted analytically.  As  a result,  the  splitter  incidence  and 
surface  velocity  distributions  are  not  accurately  predicted. 

To  obtain  a completely  satisfactory  iWDdeling  of  the  flow  in 
the  high-turning,  transonic  cascade  will  probably  require  a 
method  that  includes  both  local  viscous  effects  and  full  recogni- 
tion of  the  three-dimensionality  of  the  problem. 
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LIST  OF  SYMBOLS 


A 

AVR 

b 

C 

M 

P 

q 

(x,y) 

a 

e 

S' 

Y 

6 

0 

p 

cr 


area 

axial  velocity-density  ratio 

streamtube  height  normal  to  plane  of  calculation 

chord  length 

Mach  number 

pressure 

velocity  head 

maximum  blade  thickness 

coordinates  of  axial-tangential  coordinate  system 
angle  of  attack#  S“Y 

flow  angle  with  respect  to  axial  direction 

blade  angle  with  respect  to  axial  direction 

blade  stagger  angle  with  respect  to  axial  direction,  or 
ratio  of  specific  heats 

deviation  angle,  S,pg-S,J,g 

turning  angle, 

density 

cascade  solidity 
stream  function 


Subscripts 

1 upstrecun 

2 downstream 

/ lower  streamtube 

s static  conditions 

S splitter 

t total  or  stagnation  conditions 

TE  trailing  edge 

u upper  streamtube 

X axial  projection 
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